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ABSTRACT 


Turbulent heat fluxes in the Greenland, Norwegian and Barents Seas during March 
1988, February and March 1989, November 1991 and January and March 1992 
have been calculated with the bulk method using shipboard-based measurements 
of wind speed, air and sea surface temperatures, relative humidity and atmospheric 
pressure. The largest mean total turbulent heat flux, near 250 W/m?, was in the 
Greenland Sea in March 1989. The Norwegian Sea had mean turbulent heat fluxes 
of 130 W/m’, whereas the Barents Sea had the smallest mean turbulent heat 
fluxes. These results compared satisfactorily with climatological studies of the 
region. However, this study shows the turbulent heat fluxes to be much smaller 
than those of a recent study, especially in the northern Greenland and Barents Seas. 
Additionally, comparison of turbulent heat flux values based on 10 minute averages 
with fluxes calculated from averages of the bulk variables for an entire ship's cruise 


(10 - 22 days) shows the values to differ by only ~ 5%. 
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I. INTRODUCTION 


Determination of the turbulent (latent plus sensible) 
heat flux between the ocean and atmosphere is crucial for 
understanding several aspects of the earth's atmosphere- 
ocean system. The turbulent heat fluxes combined with the 
radiative fluxes determine the net exchange of heat across 
the air-sea interface. It is estimated that deep water 
formation in the Greenland, Norwegian and Barents Seas 
comprises about 60% of the total volume of the Arctic Ocean 
deep water (Swift, Takahashi and Livingston, 1983). 
Therefore, accurate measurements of the turbulent heat flux 
in this region are necessary to understand deep oceanic 
convection and deep water formation. 

The measurement of heat flux between the ocean and 
atmosphere is also an important factor in calculating global 
circulation of the ocean and atmosphere and therefore is 
critical to a complete understanding of the global climate. 
Heat flux is directly linked to oceanic convection and deep 
water formation in the Arctic. Formation of Arctic deep 
water causes southerly flow at depth, and the water that 
Sinks 1s replaced at the surface by the northerly flowing 
warm surface waters of the North Atlantic. This conveyor 
belt of oceanic flow is prominent in regulating the global 
climate. The warming of the ocean surface by the surface 
currents increases atmospheric instability leading to larger 
heat fluxes from the ocean to the atmosphere and increased 
cyclogenesis. The increase in storms causes more 
precipitation, lower surface salinity, increased ice 
formation and lower air temperature. As the sea ice melts 
the surface water becomes more buoyant, and the water column 
1s strongly stratified. The stratification in turn hinders 
convective overturning, slowing the formation of deep water 
and completing a negative feedback loop with a cycle of 


about 20 years (Bourke, 1994). 
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Knowledge of the amount and variability in time and 
Space of heat flux is also needed to predict the depth and 
other properties of the atmospheric boundary layer and the 
oceanic mixed layer. Turbulent heat flux and many 
properties of the atmospheric boundary layer and the oceanic 
mixed layer are used in synoptic and mesoscale atmosphere 
and ocean models. Therefore, accurate measurement of the 
fluxes offers a chance for improvement in the models. 
Additionally, an accurate account of heat flux 1s required 
to model the formation and dissipation of sea ice properly. 

Unfortunately, the Arctic has relatively few 
observations and is therefore less well understood than the 
mid-latitudes. Few measurements of sensible and latent heat 
fluxes have been made over the polar latitudes. Vowinckel 
and Taylor (1964) (hereafter VT), Bunker (1976) and Gorshkov 
(1989) have made climotological studies of the turbulent 
heat fluxes. Additionally, Hakkinen and Cavalieri (1989) 
(hereafter HC) examined turbulent heat fluxes using gridded 
Fleet Numerical Meteorology and Oceanography Center (FNMOC), 
Naval Operational Global Atmospheric Prediction System 
(NOGAPS) surface analysis for the year 1979. Their results 
gave extremely high values of total turbulent heat flux near 
the ice edge in the winter months. 

In this study, measurements of the turbulent heat 
fluxes in the Greenland, Norwegian and Barents Seas were 
based on direct measurements of bulk variables. Usually, 
computed fluxes were within the calculated error of the 
climatological studies. Additionally, comparisons were made 
with the results of calculating the sensible and latent heat 
flux based on 10 minute averaged bulk parameters versus ship 
cruise averaged bulk parameters. These results show the 
total turbulent heat fluxes are within 5% irrespective of 


averaging. 








In Chapter II, previous studies in the measurement of 
heat flux in the Greenland, Norwegian and Barents Seas are 
reviewed. In Chapter III, the data sets and methods of 
collection are described. Chapter IV contains a description 
of the bulk formula used to compute the turbulent heat 
fluxes. The results of this exercise and comparisons with 
previous studies are included in Chapter V and conclusions 


are contained in Chapter VI. 








II. PREVIOUS STUDIES 


In a review of the literature regarding heat flux, four 
sources were valuable. The first three studies are 
climatological works, whereas the final study emphasizes the 
relationship of heat fluxes to numerical ocean and 
atmosphere modeling. For each study the author’s methods, 
data sources, and assumptions made in the interpretations 


will be reviewed. 
A. VOWINCKEL AND TAYLOR 


Vowinckel and Taylor (1964) calculated long wave 
radiation, latent (evaporation) and sensible heat flux over 
the Greenland, Norwegian and Barents Seas north of 65° N. 

Their calculations are after Sverdrup (1951) 


E=k(e,-e,) V (2.1) 


LHF=E L, (232) 


where E (mm/24 hrs) is the evaporation rate; k is a 
dimensionless coefficient called the evaporation factor 
which is empirically determined from an average for the 
entire area. VT used 0.104 for k; The parameters e. (mb) 
and e, (mb) are the water vapor pressures at the sea surface 
and at the shipboard measurement height, respectively; V 
(m/s) 1s the wind velocity as measured onboard ship; LHF 
(cal / cm’ month) is the latent heat flux; and L, (cal) is 
the calculated latent heat of vaporization at that latitude. 








VT had more confidence in their measurements of 
evaporation than of sensible heat flux. Thus, they used the 


Bowen ratio to calculate sensible heat flux over the water, 


SHF Pp a 
R=—— =0.66 =~ _+ ? Ze 
LHF 1000 e,-e, ( 


SHF=R LHF (2.4) 


where R is the Bowen ratio; SHF (cal/cm’ month) is the 
sensible heat flux; p (mb) is the atmospheric pressure; T, 
(C) and T, (C) are the temperatures of the air and sea 
surface. 

The VT data for wind speed were determined based on ten 
years of synoptic data at ship M in the Norweglian-Barents 
sea. Mean wind speed was based on an evaluation of wind 
roses published by the British Meteorological Office (1959). 
Due to the paucity of data, wind Speed maps were prepared on 
a seasonal basis and monthly values were extrapolated. Air 
and sea surface temperatures were obtained from the British 
Meteorological Office (1959). Relative humidity, which is 
needed to compute e,, was taken from monthly mean maps. The 
relative humidity maps were produced based on data from ship 
M, small islands and exposed coastal stations. 

The VT flux data are provided for the Norwegian and 
Greenland Seas in five degree sections of latitude, and no 
longitudinal distinctions are made. Data were given 
monthly in units of cal/cm’? month. Also, the total heat 
fluxes here are computed by adding the sensible and latent 








fluxes. All heat fluxes shown in this thesis are in units 
of W/m’. To convert to W/m’, cal/cm’? month is multiplied by 
0.01614. Tables 2.1, 2.2 and 2.3 are VT heat flux data used 
for this study. 


a 
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Table 2.1 Latent, sensible and total turbulent heat flux 
(W/m?) over the Norwegian-Barents Sea from 80N to 75N (after 
Vowinckel and Taylor, 1964). 











Month 


January 





Table 2.2 Latent, sensible and total turbulent heat flux 
(W/m’*) over the Norwegian-Barents Sea from 75N to 70N (after 
Vowinckel and Taylor, 1964). 
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Table 2.3 Latent, sensible and total turbulent heat flux 
(W/m) over the Norwegian-Barents Sea from 70N to 65N (after 
Vowinckel and Taylor, 1964). 







B. GORSHKOV 


The World Ocean Atlas edited by Gorshkov (1983), 
contains annual and monthly maps of the various flux 
components. The atlas is in Russian with a short summary in 
English. The Russian headings, captions and other pertinent 
data were translated by Capt. Mary Lee, USAF. 

The main sources of information for compiling the maps 
were observations at 111 shore stations from 1936 to 1972. 
The methodology for flux calculations was not stated. 

The data were provided as time series depicting annual 
rates of heat flux for twelve locations in kcal/cm’ year. 
Station 10 was in the Norwegian-Barents Sea, station 11 was 
in the Greenland Sea and station 12 in the southern 
Norwegian Sea (refer to Figure 2.1 for locations; all 
figures are at the end of chapters). To convert to W/m, 
kcal/cm* month was multiplied by 16.14. Tables 2.4, 2.5 and 


2.6 contain heat flux data germane to this study. 











Table 2.4 Latent, sensible and total turbulent heat flux 
(W/m?) at station 10 near 72N/020E over the Norwegian- 
Barents Sea (after Gorshkov, 1983). 










LHF SHF THF 


Table 2.5 Latent, sensible and total turbulent heat flux 
(W/m?) at station 11 near 78N/O010E over the Greenland Sea 
(after Gorshkov, 1983). 










Table 2.6 Latent, sensible and total turbulent heat flux 
(W/m?) at station 12 near 64N/005E in the Norwegian Sea 
(after Gorshkov, 1983). 








C. BUNKER 


Bunker (1976) computed energy flux over much of the 
world, including one station in the northern Norwegian and 
Barents Seas that was useful for this thesis (see Figure 
2A): 

Bunker calculated heat flux based on the bulk 
aerodynamic method. His method used exchange coefficients 
that vary with wind speed and stability. His equations were 


in the form 


LHF=p Cyl, (O5-Q19) Uy (2-5) 


SHF=9CyCp(Ts-Tig) Uy (2.6) 


where p is atmospheric density; C, and C, are the variable 
exchange coefficients for water vapor and sensible heat; Ce 
is the specific heat constant; OQ, and Q,, are the averages of 
the mixing ratio of air in contact with the surface and at 
10 m; U,, is the average wind speed at 10 iy and. 72 ‘and. 25, 
are average temperatures of the sea surface and air at 10 m. 

Bunkers' values for the exchange coefficients were 
obtained from tables calculated from experimental data. The 
values of the coefficient for water vapor have been applied 
to the computations of sensible heat, making the 
coefficients equal. 

Bunkers' calculations were based on a National Climatic 
Center data set, which was collected for the years 1941- 
Lg72 . Averages for each month of each year for an entire 
Marsden Square were formed for fluxes and basic 


meteorological variables. Monthly and annual averages for 


10 








the 32-year period were formed for subdivisions of the 
Marsden Squares. A Norwegian-Barents Sea area was selected 
at 71 N, 17 E (refer to Figure 2.1 for location). MThis data 
set was produced based on 1859 observations over the 32-year 
period. Table 2.7 portrays the heat flux for this station. 


Month LHF SHF THE 





Table 2.7 Latent, sensible and total turbulent heat flux 
(W/m*) over the Norwegian Sea at 71N/017E (after Bunker, 
1976). 


D. HAKKINEN AND CAVALIERI 


Hakkinen and Cavalieri (1989) estimated oceanic surface 
heat fluxes in the Norwegian, Greenland and Barents Seas 
uSing gridded Fleet Numerical Meteorology and Oceanography 
Center (FNMOC) surface analysis of winds and temperatures 
for one year. They used the bulk method (Equations 2.2 and 
2.3) with the heat (C,), and moisture (C,) transfer 
coefficients fixed at a value 1.5 x 10°? neglecting any 
Stability dependence. Density was also fixed at 1.3 kg / 
m°. 

The data for their calculations are based on archived 
Navy Operational Global Atmospheric Prediction System 
(NOGAPS) analyses for the period December 1, 1978 through 


November 30, 1979. The gridded data has a resolution of 250 


ad 








km. The NOGAPS model data includes sea level pressure, air 
temperature, vapor pressure, winds and SST. All parameters 
except SST are given twice daily. SST is given once daily 
and interpolated to 12 hour intervals for flux computations. 
Monthly mean contour plots of total turbulent (sensible 
and latent) heat fluxes were provided in W/m’. Tables 2.8 
and 2.9 contain total turbulent heat fluxes (LHF + SHF) as 
interpreted from contoured maps for regions and times 


corresponding with data analyzed in this thesis. 







Month 


THE 
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Table 2.8 Total turbulent heat flux (W/m?) over the 
Greenland Sea 77.5N/005E (after Hakkinen and Cavalieri, 
1989). 


a aa 


Table 2.9 Total turbulent heat flux (W/m?) over the 
Norwegian Sea 68N/005E (after Hakkinen and Cavalieri, 1989). 
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Figure 2.1 Location map for climatological stations. 


Stations 10, 11 and 12 are data from Gorshkov, and station B 
1s data from Bunker. 
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III. DATA 


A. COLLECTION OF DATA 


This study utilizes data collected by Professor Kenneth 
Davidson and associates at the Naval Postgraduate School 
(NPS) in the Greenland, Norwegian and Barents Seas during 
several recent ship cruises (Figure 3.1). A Coastal Climate 
WeatherPak meteorological station onboard the Research 
Vessel (R\V) Haakon Mosby sampled air temperature, wind 
Speed and direction relative to the ship's heading every 
second and relative humidity and atmospheric pressure every 
12 seconds. The measurement platform was at 15 m above the 
sea surface. In addition, the ship’s speed and heading were 
recorded, enabling the calculation of true wind speed and 
direction. Sea Surface Temperature (SST) was also measured 
by the ship. All observations were averaged over 10 minute 
intervals. Details of the various ship cruises are 
described in later sections. 

The temperature data were estimated to be accurate to 
within 1.0 C and the relative humidity data to within 53%. 
Wind speeds were estimated to be accurate to within 0.3 m/s. 
When winds were from the stern, errors are probably greater 
due to the sheltering effect caused by the ship’s 
superstructure. The measurements of atmospheric pressure 
are accurate to within 2 millibars. 

At times, data were missed or not recorded; during 
these cases the trends were assumed to be linear and 
intermediate data points were interpolated. 

A common bias involving data collected by ships is 
Storm avoidance. That is, if a ship's crew expects a 
location to be stormy they will typically avoid that area. 
This 1s not so with the cases studied here. In general, the 
locations and ship tracks were planned in advance and were 


not significantly deviated from. The only storm avoidance 


i> 








bias occurred during the later part of SIZEX in January 
1992, when an intense cyclone disabled the weather 
collection unit with winds greater than 28 m/s. This 
prevented about 24 hours of data collection during a high 
wind period. This was estimated not to Significantly affect 


the overall heat flux averages. 


NORCSEX | CEAREX NORCSEX SIZEX 
Norw Barents Norw Barents 


Mar 88 Feb 89 Nov 91 





Table 3.1 Mean and standard deviations of input data for 


indicated exercises. 


B. NORCSEX 88 


The Norwegian Continental Shelf Exercise (NORCSEX) 88 
was a Buropean Research Satellite (ERS-1) pre-launch 
calibration and validation experiment. The purpose of the 
ship cruise was to validate data measured by aircraft flown 
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with ERS-1 sensors. NORCSEX 88 was conducted in the 
Norwegian Sea and centered around the Norwegian Continental 
Shelf at 64°30’N and 009°E during March of 1988. From 07 to 
18 March, the R/V Haakon Mosby collected meteorological data 
while transiting along the west coast of Norway (see Figure 
3.1 for track). The time series of these data are shown in 
Figure 3.2. 

NORCSEX 88 displayed the lowest wind speeds of all the 
cases; during the exercise moderate winds prevailed with a 
mean of only 7.0 m/s. However, gale to storm force 
conditions occurred on 11 March (Julian day 71); as the wind 
veered from the south to the north, the wind speed quickly 
increased to more than 22 m/s. Along with the wind change 
came a decrease in air temperature and an increase in 
relative humidity. Throughout the measurement period, the 
alr temperature and sea surface temperature were fairly 
constant. SST had a range of 2 C to 8 C and a standard 
deviation of 1.1. Air temperature had a range of -5 C to 5 
C and the smallest standard deviation of all the exercises 
Aw wa: (Ce 


C. CEAREX 


The Coordinated Eastern Arctic Exercise (CEAREX) was a 
large international interdisciplinary project involving 
several ships, aircraft and ice camps. CEAREX was conducted 
in the vicinity of the Svalbard Islands from September 1988 
to May 1989. For this study, continuous surface 
measurements collected by the R/V Haakon Mosby from 25 Feb 
to 23 Mar 1989 were used. From 25 February to 09 March 1989 
the ship was in the northern Norwegian and Barents Seas 
(Figure 3.1). Since these tracks correspond more closely to 
land and climatological stations in the Barents Sea, this 
region will be called only the Barents Sea to prevent 


confusion with locations in the southern Norwegian Sea. 


di, 








After 10 March, the ship transited on a westerly course to a 
location west of Svalbard and conducted north-south tracks 
along the ice edge. Since the measurements in the Barents 
Sea differed markedly from those in the Greenland Sea, the 
measurements have been divided into two separate data sets. 
The time series of data in the Barents Sea is shown in 
Figure 3.3, and the time series for data in the Greenland 
Sea 1s shown in Figure 3.4. 

There were several cyclone passages while the ship was 
located in the Barents Seas. The temperature tended to be 
near freezing for relatively long periods of time with brief 
surges of extremely cold air after cyclone passages. The 
average wind was strong at 10.7 m/s and primarily from the 
north. Sea surface temperatures were not available for 
CEAREX. Mean SSTs of 1.8 C and 1.4 C were chosen for CEAREX 
data in the Barents and Greenland Seas, respectively. These 
values were based upon averages from SIZEX data along 
Similar tracks and climotological data. 

In the Greenland Sea, the air temperatures were much 
colder than in the Norwegian Sea, with mean temperatures of 
-9 C compared with -1 C. After passage of a cyclone on 
about 11 March (Julian day 69) the pressure field was 
relatively flat with colder temperatures and lower relative 


humidity. 
D. NORCSEX 91 


The Norwegian Continental Shelf Exercise (NORCSEX) 91 
was a follow-on to NORCSEX 88. Its purpose was a post 
daunch validation of ERS-1 sensors. NORCSEX 91 was 
conducted in the Norwegian Sea from 08 Nov to 29 Nov 1991. 
During this period meteorological variables were collected 
onboard the R/V Haakon Mosby. The time series of data are 
shown in Figure 3.5. 

NORCSEX 91 data contained the warmest air and sea 
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surface temperatures of all the exercises with means of 5.8 
C and 8.3 C, respectively. The temperatures were also less 
variable than those of other exercises with standard 
deviations of 2.2 C and 0.5 C. During the period there were 
several strong cyclones that passed over the ship. These 
cyclones exhibited higher winds but showed little 


temperature change. 


E. SIZEX 


The Seasonal Ice Zone Experiment (SIZEX) was conducted 
as a part of an ERS-1 ice Synthetic Aperture RADAR (SAR) 
validation. SIZEX was held in the Greenland Sea during 
1992. The experiment has been divided into two periods, one 
in January and one in March. Time series of data are shown 
in figures 3.6 and 3.7. 

In January, the data were influenced by three prominent 
cyclone passages. The cyclones were accompanied by strong 
winds from the south to southeast. The final cyclone, 
occurring on about 13 January (Julian day 13), was storm 
force with wind speeds approaching 30 m/s. The storm was 
strong enough to knock out the weather collection station. 

In February and March, there were no intense cyclones 
but periods with significant winds between 15-20 m/s. 
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Figure 3.1 Ship tracks of NPS exercises. 
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Figure 3.2 Time series of wind Speed, wind GlLrection: 
Surface Cemperature (dotted) and air tCemperature 


sea 
relative Humidity and Pressure during NORCSEX88, 


(solid), 
March 1988. 
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Figure 3.3 Time series of wind speed, wind direction, sea 
surface temperature (dotted) and air temperature (SOl1d) 
relative humidity and pressure during CEAREX in the Barents 


sea, February-March 1989. 
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Figure 3.4 Time series of wind speed, wind direction, sea 
Surface temperature (dotted) and air temperature (solid), 
relative humidity and pressure during CEAREX in the 
Greenland Sea, March 1989. 


o a 
E 90 | 
210 1? \ 
2G 
315 320 305 330 335 
3 eas . ‘. 
v J ; 
ey | aa 
: 315 390 325 330 335 





315 320 325 336 335 
100 
E sol OY 
0 Sl Spm i ee a 
315 320 325 330 335 
1020 
< 1000 wre 
a 9S | 
oo 315 320 325 330 335 
julian date 


Figure 3.5 Time series of wind speed, wind direction, sea 
surface temperature (dotted) and air temperature (solid), 
relative humidity and pressure during NORCSEX91 in the 
Norwegian Sea, November 1991. 
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Figure 3.6 Time series of wind speed, wind direction, sea 
Surface temperature (dotted) and air temperature (solid), 
relative humidity and pressure during SIZEX in the Greenland 
sea, January 1992. 
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Figure 3.7 Time series of wind speed, wind direction, sea 
Surface temperature (dotted) and air cemperature (solid), 
relative humidity and pressure during SIZEX in the Barents 
sea, March 1992, 
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IV. HEAT FLUX COMPUTATIONS USING THE BULK METHOD 
A. THEORY 


The primary goal of this thesis is to quantify the 
turbulent heat flux between the ocean and atmosphere in the 
Greenland, Norwegian and Barents Sea regions. This goal is 
accomplished through use of the bulk aerodynamic method 
which is based on Monin-Obukhov (M-O) similarity theory. 

The bulk method has proven to be an accurate method for 
computing heat fluxes between the ocean and atmosphere 
(Blanc, 1985). This method is the most practical for 
shipboard-based measurements, requiring only routine 
meteorological measurements such as wind speed, air 
temperature, sea surface temperature, relative humidity and 
atmospheric pressure. The bulk method allows calculations 
to be made based on measurements made from ships at sea, 
which provide a substantial portion of the data available 
over the oceans. 

Monin-Obukhov Similarity theory assumes the gradients 
of the meteorological variables depend on surface fluxes, 
height and stability only. For instance, the wind is 
modeled using the log wind profile modified for stability 
effects. Therefore, if the wind speed is known at any 
height within the surface layer, one can estimate the wind 
speed at any other height within the surface layer (Equation 
4.1). Likewise, the profiles of temperature and moisture 
can be used to estimate sensible and latent heat fluxes. 
The profiles for temperature and humidity are known and can 


be found by Equation 4.2, 
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Where u(m/s), T (C) and z(m) are wind speed, potential 
temperature and observation height; u. (m/s) is a velocity 
scale called the friction velocity, T, (C) is a surface 
layer temperature scale; k is Von Karman's constant; and @Q, 
and ®, are dimensionless wind shear and temperature 


functions of stability (z/L). 
B. BULK FORMULA 


The turbulent heat fluxes at the surface can be 


expressed in the form 


SHF =p c,w'T!= -p Chu, T, (4.3) 


LHF=p L,w'q’=-pL,u, q, (4.4) 


where SHF (W/m*’ ) and LHF (W/m?) are sensible and latent 
heat flux, with the sign convention that a positive value 
indicates a heat flux from the ocean to the atmosphere and a 
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negative indicates heat flux from the air to the ocean; The 
variable p (kg/m’)is atmospheric density, based on surface 

temperature and pressure and calculated using the ideal gas 
law; C, (1004 J/kg K) is the specific heat of dry air; L 
(250000 J/kg)is the latent heat of vaporization at 0°C. 


Vv 


The parameter T. (C) is a temperature scale and q. (g 
water/g air) is a humidity scale. The scaling paramters are 
calculated using the following equations 


uk 


in(—)-¥, 
Zo 


u, (4.5) 


ae k 
In(—)- ¥, 


Zot 


T, = (4.6) 


g, =! ~Fs) K (4.7) 


1n(—+)-¥, 
Zog 


where Ti. (C) and T, (C) are potential temperature at 10 m 
and at the surface; q,, (g/kg)and q, (g/kg) are the mixing 
ratios at 10 m and at the surface; z, is roughness length, 
which is the sum of Charnock’s roughness length and the 
roughness length for a smooth surface and is defined by 


Equation 4.8; Zz, and Zo, are the thermal and moisture 
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roughness lengths based on a constant heat exchange 


coefficient (C,, and C,, ) and defined by Equations 4.9 and 
4.10; 


Z)=(0.011u,s/9.8) +( (.11) (1.4x107) /u,) (4.8) 


-k? 
Zop=10 Exp—__—_—__—_ (4.9) 
ue agro 1n(z/2Z,) 

~[p2 
Zog = 10exp . 4.10 


Cop in(2z/Z,) 


¥, and ¥, are the integrated stability functions, for 


unstable stratification and defined by Equations 4.11 and 
4.12. 


F,=2 In[(1+x) /2] +1n[(1+x?) /2] -2 arctan(x) +x/2 (4.11) 


¥.=2 In[(1+x?) /2] (4.12) 


X=(1-16 2/1) -35 (4.13) 
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and tor stable stratification 


Piz-5 z/L (4.14) 


P.=-5 z2/L (4.15) 


L is the Monin-Obukhov stability length defined by 
Equation 4.16, and z/L is the stability; Ty. 1S the buoyancy 
flux scaling parameter determined by Equation 4.17, 

u,? T. 


-_* "Vs .16 
vis kot (4 ) 


(T,-T..) K 


“"Tnlz/z,) -8,  '4:17) 


where T, (C) and T,, (C) are the virtual potential 
temperature at the measurement height and the surface. 

The constants in the Charnock relationship (Equation 
4.8) and the integrated stability functions are from Smith 
(1988) . 

The advantage of using the bulk method is that the only 
parameters required are the mean wind speed, the mean air 
temperature and the mean humidity at one level each in the 
surface layer, and the mean surface temperature. High 


frequency or multi-level measurements are not required and 
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the bulk method is less affected by ship effects than are 
other flux calculation methods. 

For the unstable cases (cold air over warm water) and 
stable cases (warm air over cold water) W is a function of 
the stability parameter (z/L). The surface layer over open 
water in the Arctic is typically unstable. Therefore, z/L 
is negative, and the sensible and latent heat fluxes are 
enhanced by buoyancy effects. However, if the surface layer 
1s stable, which is rare but does occur, then z/L is 
positive. Thus, VP, is negative and u., T. and q. would all 
be smaller in magnitude; this implies the sensible and 


latent heat fluxes would also be smaller in magnitude. 
Cy SOLVING FOR HEAT FLUX BY ITERATION 


Determining us, T., and q. requires knowledge of L, 
which in turn is a function of u., T. and q.. Therefore, an 
iterative solution method is needed. 

The first step in the iterative process is to make an 
initial guess at L, Zo, and us. Then Equations 4.5, 4.6, 
4.7, 4.8 and 4.16 are repeatedly solved for u,, ZeJawe 
until the solution for u. converges. However, if the wind 
speed 1s zero or the atmosphere is too stable, the solution 
will not converge. In this case the atmosphere is non- 
turbulent and u., T., gq. and z/L are set equal to zero. The 
Matlab program used in calculating heat fluxes can be found 


in the Appendix. 
D. LIMITATIONS OF THE BULK METHOD 


The scheme shown above is one of many bulk transfer 
coefficient schemes. Each of the different schemes contains 


different coefficients of drag, heat and humidity and varies 
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in degree of complexity. Some include stability dependence 
while others ignore stability or assume a constant stability 
effect. Additionally, each method has a different range of 
wind speeds and air sea temperature differences under which 
the method is considered valid. Smith (1988) listed several 


limitations to his technique, which are summarized below. 


1. High Wind Speed 

There are few reliable measurements of wind stress and 
heat flux above 26 m/s. Thus, application of the method to 
wind speeds above 26 m/s would be an extrapolation of the 
lower wind speed data. Fortunately, winds above the 26 m/s 
threshold are absent in all but the SIZEX data set and are 


rare there. 


2. Extremely Stable Stratification 

If the atmosphere is too stable then the iterative 
method will fail to produce a convergent solution. Here the 
program will give a value of zero for u*, Te>. G an 27 lin 
This implies the fluxes are zero and actual values will be 
very close to zero. These highly stable conditions would 
occur if the air temperature was very warm and the sea 


surface very cold, but such occurrences are rare. 


3. Extremely Unstable Stratification 

If the wind speed were close to zero then the iterative 
method would fail to produce a convergent solution. In this 
case the program will give a value of zero for u*, which 
implies the fluxes are zero. In fact, there may be fluxes 
due to convection. Free convective scaling was not included 
in the bulk formula, but this was judged to not degrade the 
results since calm conditions and non-convergent solutions 


are very rare or nonexistent in the data. 


a 
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V. RESULTS AND COMPARISONS 


A. RESULTS OF HEAT FLUX COMPUTATIONS 


Turbulent heat fluxes were computed based on the data 
presented in Chapter III and the method discussed in Chapter 
IV. Overall, the average measured turbulent heat fluxes 
(THF) in this thesis (Table 5.1) compared favorably with 
those given by the climatological studies mentioned 


previously. 






Exercise 
Date 


Region 











CEAREX 
Mar 89 
Greenland Sea 








SIZEX 
Jan 92 
Greenland Sea 








NORCSEX 
Mar 88 
Norwegian Sea 











NORCSEX 
Nov 91 
Norwegian Sea 











CEAREX 
Feb 89 
Barents Sea 








SIZEX 
Mar 92 
Barents Sea 






Table 5.1 Mean value of latent, sensible and total 
turbulent heat flux (W/m?). 


ie. 








Heat flux computations for all the NPS ship cruises 
occurred during the Arctic winter (November - March). 
Analysis of the climatological data of Bunker (1976) and 
Gorshkov (1988), showed peak heat loss (positive or upward 
heat flux) for this region in December and January with only 
slightly less heat loss in the preceding and following 
months. Therefore, the time difference in the observations 
does not significantly hinder comparisons between exercises 
that occurred in different winter months, and the 
differences in heat flux values displayed in Table 5.1 are 
primarily due to the differenecs in the ship's location or 
differences between years. 

According to the data, the largest mean THF, near 250 
W/m’, was seen in the Greenland Sea. The Norwegian Sea had 
a mean THF of 130 W/m’, whereas the Barents Sea had the 


lowest mean THF. 
B. COMPARISON WITH PREVIOUS STUDIES 


In this section, the calculations of heat flux based on 
the NPS data are compared with previously published values. 
Accepted faults with comparisons include: (1) Slightly 
differing spatial regimes: the mean value of the total 
turbulent heat flux measurement is taken along a ship track, 
whereas in the climatological studies locations are fixed or 
averaged over large areas and HC used grid point data. (2) 
Different time scales: record length was two to three weeks, 
whereas the climatological studies have tens of years of 
data, and HC data comes entirely from one year (1979). 
However, the bulk method used here, is similar to Bunkers! 
method, and is considered more advanced than some earlier 


methods. For example, the method used in this thesis 
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includes flux computations that account for variable 

stability (z/L), roughness (Zo) and atmospheric density (p). 
In the following subsections, the results of each 

exercise's heat flux computations are reviewed and compared, 


by region, to the results of the earlier studies. 


1. Greenland Sea 

During SIZEX, the mean LHF was 111 W/m? and the mean 
SHF was 105 W/m*. This exercise had some of the largest 10 
minute averaged heat fluxes of all the cases (Figure 5.1). 
These large heat fluxes are a result of the stronger winds 
and strong to moderate air sea temperature differences 
(ASTD) encountered in the Greenland Sea during SIZEX. 
Histograms of THF for all exercises were plotted with 50 
equally spaced bins between the maximum and minimum values 
of THF for that exercise. The histogram for this case 
(Figure 5.2) showed a maximum occurrence near the median of 
210 W/m’ and a long tail to the right (above 400 W/m?) 
corresponding to a high wind event on Julian day 14. A 
false local maximum occurring at 100 W/m? is the result of a 
data gap of about six hours on Julian day 13. 

During CEAREX, the heat fluxes (Figure 5.3) were also 
relatively large compared with other NPS exercises. The 
mean LHF was 109 W/m* and the mean SHF was 140 W/m’. As in 
SIZEX, the large heat fluxes can be attributed to large 
ASTD, large moisture gradients and moderate winds. The 
large ASTD values are a result of the cold air advection off 
the ice (northerly winds). Mean winds (9.1 m/s) for the 
Greenland Sea region during CEAREX are similar to other 
exercises and had a small standard deviation (2.7 m/s) 
compared with other exercises. Nevertheless, the heat 
fluxes are larger in the Greenland Sea primarily due to the 


colder air over the Greenland Sea. The histogram of the 
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heat flux (Figure 5.4) has at least three peaks. The 
different peaks are the result of changes in synoptic 
conditions. The largest number of occurrences is at 330 
W/m’ and is related to the cold air outbreak on Julian days 
7i through 75. During this period the highest heat fluxes 
were caused by cold air advection, and occurred after the 
winds backed to the north and result in a highly unstable 
Surface layer (very cold air over warmer water). The second 
peak is near 200 W/m* and is the result of smaller winds and 
ASTD after Julian day 75° The third local maximum at 60 
W/m*, occurred during a pressure minimum, and was the result 
of the southerly winds that gave an unstable surface layer 
(cold air over warmer water). 

The VT heat flux data compared well with the January 
SIZEX data, but were 20% larger than the March CFAREX 
observations of heat flux (Tables 5.2 and 5.3). The 
Gorshkov data showed mixed results. Their values for heat 
flux during SIZEX were within 10% and during CEAREX were 702 
less than computed for the NPS data set. The HC heat fluxes 
were two to three times larger than any other computed heat 
flux during January, and were 20% larger than the NPS values 


during CEAREX, which was in March. 


eis 








Table 5.2 Latent, sensible and total turbulent heat flux 


HC 





(W/m*) over the Greenland Sea in January. 


CEAREX GS 109 





Table 5.3 Latent, sensible and total turbulent heat flux 


(W/m*) over the Greenland Sea in March. 


2. Norwegian Sea 

NORCSEX 88 was conducted in March 1988, off the 
southwest coast of Norway. Heat fluxes for this exercise 
were Similar to previous results for the area (Table 5.4). 
The mean LHF was 87 W/m* and the mean SHF was 44 W/m. 

The heat flux (Figure 5.5) exhibits a synoptic 
influence, which should be expected in a stormy region of 
the Norwegian Sea. For example, before the passage of a 
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series of surface cyclones on Julian day 71 the THF is small 
( 50 W/m* ) due to low winds and a small ASTD. Afterwards, 
the wind speed and ASTD increased causing the THF to climb 
above 200 W/m’ and then decreased to around the mean of 130 
W/m’ after the winds died down. Also, the SHF time series 
1s strikingly similar to the time series of wind speed (as 
seen in Chapter III, Figure 3.1). The histogram of heat 
flux (Figure 5.6) is nearly symmetric with a sharp peak near 
the median of 120 W/m’. The tail of fluxes above 250 W/m’ 
represents a mesoscale gale force storm observed for several 
hours on Julian day 71. 

NORCSEX 91 was conducted in the southern Norwegian Sea 
in November 1991. The THF values for this exercise (Figure 
5.7) were similiar to those of NORCSEX 88. During NORCSEX 
91, the SHF was low with a mean 30 W/m’. The LHF was 
larger, with a mean value of 96 W/m’®. NORCSEX 91 
experienced the passage of several cyclones with associated 
wind variations. The heat fluxes modulated with these wind 
changes. Also, apparent is the effect of the ASTD. During 
the earlier periods with larger ASTD the fluxes are larger 
but after Julian day 328 when the ASTD became smaller the 
flux decreased. The histogram of THF (Figure 5.8) has three 
modes. The first mode, near 0 W/m’, was a result of the 
small or zero ASTD that occurred from Julian day 353-00 335. 
Modes two and three occurring at 90 W/m? and 170 W/m are 
the result of a bimodal distribution about the median of 125 
W/m: 

The differing data sets compare well within this region 
as all the data sets are within 20% - 30% of the mean 
observed THF (Tables 5.4 and 5.5). 
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Table 5.4 iaweue: sensible and total turbulent heat flux 


(W/m) over the Norwegian Sea in March. 


NORCSEX 91 















a a 
a CC CC 


Table 5.5 Latent, sensible and total turbulent heat flux 





(W/m*) over the Norwegian Sea in November. 


3. Barents Sea 

A segment of CEAREX was conducted in the Barents Sea in 
February and early-March 1989. During CEAREX in the Barents 
Sea the mean LHF was 73 W/m’? and the mean SHF was 45 W/m’. 
several cyclones seen in the data are apparent in the time 
series of heat flux as peaks at intervals of two to three 
days (Figure 5.9). The histogram of THF (Figure 5.10) is 
Skewed to the right and has some high peaks at 80 W/m. 
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Also noticeable are the negative fluxes (heat flux into the 
ocean) resulting from the surge of warm southerly air and 
stable conditions around Julian day 66. 

SIZEX was conducted in the Barents Sea in March. The 
THF values for March in the Barents Sea were small, near 60 
W/m’. During SIZEX, the LHF was 44 W/m’ and the SHF was 20 
W/m’. The heat fluxes are smaller here due to smaller ASTD 
(Figure 5.11). The histogram of THF (Figure 5.12) is 
symmetric with a median of near 60 W/m’. 

The heat fluxes of VT and Bunker are up to two times 
greater than observed by NPS ship data in the Barents Sea 
(Tables 5.6 and 5.7). However, the Gorshkov and HC heat 
fluxes are three to four times larger than the measured 
values. The values of heat flux in the NPS data are 
strongly influenced by the lack of SST data and the 
assumption of a constant 1.8 C SST for CEAREX. For example, 
1f the average values for all other directly measured 
quantities (wind speed, air temp, rh, p) were unchanged and 
the average SST was assumed to be 0 C vice 1.8 C, the THF 
would become 60 W/m’, nearly half. If the average SST were 
> C, the THF would become 230 W/m? or nearly double the 
estimated value. Normally SST would not be so important. 
However in this instance it is critical because the assumed 
SST (1.8 C) and the actual mean air temps (-1.4 C) are 


Similar. 
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a a Ce 
a LC 


Table 5.6 Latent, sensible and total turbulent heat flux 


(W/m’) over the Barents Sea in February. 


Gz CC Cc 
a Ce 
ae Cc 

Table 5.7 Latent, sensible and total turbulent heat flux 


(W/m*) over the Barents Sea in March. 
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C. LIMITATIONS OF PREVIOUS STUDIES 


1. Time Averaging Methods 

Hanawa and Toba (1987) (hereafter HT) made comparisons 
of data averaged over various time periods from one day to a 
month using the sampling method (SM) and the scalar 
averaging method (SAM). The sampling method (SM) is a 
Simple average of fluxes computed from hourly (or every 
observation) values of the bulk meteorological parameters 
(wind, temperature and humidity). The scalar averaging 
method computes the average flux from average values of the 
bulk data, thus ignoring correlations between variables. 
Both methods were compared with data taken every three hours 
from a mid-latitude Oceanographic Weather Station (OWS-T) 
(29N 135E) from June 1950 to November 1953. HT showed that 
the SHF was almost the same regardless of method and 
averaging time. Also, LHF values computed by the SAM are 
about 105% of those by the SM for averaging times from three 
days to one month. 

Ledvina et al. (1993) made comparisons similar to HT, 
using averaging times from 2-72 hours. Their ratios vary 
little for times exceeding 36 hours, thereby eliminating the 
need to extend the averaging period beyond 72 hours. Their 
data was in the equatorial Pacific near O000N 145E from 17 
February to 10 March 1990. Like HT, they showed the SAM 
overestimated LHF with increasing averaging period. Yet 
with results drastically different from HT, they showed the 
SAM to underestimate SHF by 40%. Much of this difference 
may be related to convection in the tropics, thereby 
demonstrating the importance of regional influences on heat 
fluxes. 

The purpose of this subsection is to extend the work 
completed by HT and Ledvina et al. (1993) to the polar 
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latitudes with recent data, thus providing a basis for 
determining the accuracy of heat fluxes computed via various 
averaging techniques in the polar seas. Ideally, time 
averages are desired that are long enough to include all low 
frequency effects (synoptic) and sampled fast enough to 
capture all high frequency effects (Kaimal and Finnigan, 
1994). 

Esbensen and Reynolds (1981) define the sampling method 
for calculating SHF and LHF in the following equations 


SHF = pc,C,0AT (5.1) 


LHF = pL, C, U Ag (5.2) 


where averaging is indicated by the overbar. The key 
feature of this method is that heat fluxes are calculated 
for every measurement. The measurement interval was 10 
minutes in the NPS data set, generally 3 hours in Bunker's 
data to 12 hours in the HC model. The fluxes are then 
summed and then averaged. 

The scalar averaging method (classical method) for 
computing heat flux defines the flux in the following forms 


SHF = pC,C, UAT (5.3) 


LHF =p L,C,UAq (5.4) 


where the average transfer coefficients of heat and humidity 
are computed empirically from the sample means of 
temperature difference, humidity difference and wind speed. 
The difference of the classical method (SAM) from the 
sampling method is the measurements are averaged, usually 
over the entire month and then the heat fluxes are computed 
based on the averages of the measurements. 

Lumley and Panofsky (1964) suggest that a reasonable 
averaging time for average mean horizontal winds (and 
temperature) would be 30 minutes, corresponding to the 
passage of two to three large convective cells through the 
depth of the convective boundary layer. However, if 
averaging periods are increased to about one hour, 
nonstationarity in the form of diurnal variations will 
occur. With averaging periods of one hour or greater, 
errors. Of. 5% and “L0=50%- for ‘TT. and.w result. 

Bunker (1976) and HC used the sampling method when 
calculating their heat fluxes; VT and Gorshkov (1988) used 
the classical approach or scalar averaging method. Finally, 
data from this study are used to compare the sampling method 
with the classical method in the Greenland, Norwegian and 
Barents Seas. 

The approach of Bunker (1976) was to include every ship 
or buoy report and compute heat fluxes and then average the 
fluxes over the month. This method was better than the 
Classical approach in that it would not result in errors due 
to diurnal fluctuations. The classical approach used by 
Gorshkov (1988) and VT also has it advantages. Since the 
covariance and correlation of UAT and UAg are small at 
shorter time intervals and increase for longer periods, the 
classical method may provide a better heat flux estimate 
when the data set is small (Ebsensen and Reynolds, 1981). 


The smaller covariance at shorter time periods probably 
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comes from the fact that stronger winds cause more mixing 
and thus decrease AT and Aq. 

Results from NPS cruise data show the SAM (using 
average wind speeds, temperatures, relative humidities and 
Pressures for the entire period of the exercise) gave heat 
flux results that are ~5% (6% LHF and 33 SHF) larger than 


when using the sampling method (Table 5.8). 






Exercise 






Region 






CEAREX 
Greenland Sea 







SIZEX 


Greenland Sea 








NORCSEX 88 


Norwegian Sea 







NORCSEX 91 


Norwegian Sea 






CEAREX 
Barents Sea 








SIZEX 


Barents Sea 






Table 5.8. Ratio of heat flux (W/m?) calculated using time 


averages of meteorological variables versus uSing the 


sampling method. 
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2. Stability Dependence 

If stability is neglected, then one must assume neutral 
or near neutral conditions and buoyant convection would be 
negligible. However, neutral stability is rare and unstable 
conditions prevailed during the NPS exercises. The effect 
of stability 1s applied through use of the stability 
functions (W, and YW.) defined in Equations 4.11 and 4.12 for 
unstable stratification and Equations 4.14 and 4.15 for 
stable stratification. For unstable stratifications 
(negative z/L) both stability functions (¥) are complex 
functions of stability, and their values will be relatively 
small and negative thus yielding larger LHF and SHF than 
under neutral conditions. For stable stratification 
(positive z/L), Wis a linear function of stability (W = - 
9Z/L) and will always yield smaller LHF and SHF than neutral 
conditions. 

The method used by Bunker accounts for stability. HC 
however, do not explicitly include stability in their 


Calculations. 


3. Transfer Exchange Coefficients 

The average C, and C, for the NPS Barents Sea exercises 
were 1.35x10°* and 1.13x10°, respectively. If Bunker's 
nomogram for transfer coefficients were used instead, that 
value would be 1.58x10°°. This 17% and 40% increase in the 
exchange coefficients would lead to a corresponding increase 
in LHF and SHF. Bunker's C,, and C,, values for near-neutral 
conditions were strongly dependent upon wind speed and AT. 
This dependence on wind speed is however unwarranted as 
there does not seem to be conclusive experimental evidence 
for variation of C, and C, with wind speed (Smith, 1988). 

HC heat flux errors would result from the use of 


constant exchange coefficients, vice coefficients that vary 
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with wind speed, temperature gradient, humidity gradient and 
stability. ‘HC £ixed.C. and -C. 6b 1.510 for all. regions. 

In contrast, the coefficients used here varied by region 
with C, having a range of 1.0x107 to 1.25x10° with an 
average of 1.13x10° and C, having a range of 1.2x10%? to 
1.5x10°* with an average of 1.35x10°. I£ calculating fluxes 
for our exercises, this variation would cause HC to 
overestimate SHF nearly 45% in the Barents Sea and 27% in 
the Norwegian and Greenland Seas and to overestimate LHF 253% 
in the Barents Sea and 9% in the Norwegian and Greenland 


Seas. 


4. Other Factors 

As discussed in Chapter II, the VT data are divided 
into five degree latitudinal sections and no distinctions 
are made for longitude. Based on this I have labeled the VT 
75N - 80N region as the Greenland Sea, their 75N - 70N as 
the Barents Sea and their 70N - 65N as the Norwegian Sea. 
These distinctions are not absolutely accurate ina 
geographic sense, but they are fitting with the 
climatological stations of Gorshkov and Bunker and with our 
ship cruise data. 

Haikkinen and Cavalieri used gridded data for the 
calculation of heat flux. Thus, comparison of locations 
should be easily accomplished. However, their data was for 
only one year (1979), which will likely cause errors when 
comparing with different years. Other errors MiGue: OCCur. by 
use of model data vice observed. If the model was biased 
toward higher winds or colder air temperatures then larger 
heat fluxes would be shown. Additionally, the coarse model 
grid (250km) does not have the resolution to accurately 
portray the large horizontal gradients of temperature and 


wind speed that occurs at the ice edge (Personal 
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communication Guest, 1995). This may in part account for 


their large flux gradients near the ice edge. 
D. POWER SPECTRAL DENSITY OF HEAT FLUXES 


Power spectral densities (psd) of heat flux for all 
exercises were also examined. The power spectral density 
function estimates the signal of heat flux as a function of 
frequency. Windowing was performed using a Hanning window. 
The sampling frequency that varied for each exercise based 
on the record length was used for scaling plots. One half 
of the sampling frequency (1 / 600 sec ) would be the 
highest resolvable frequency (Nyquist frequency). The 
lowest resolvable frequency would be 1 / the total time of 
each exercise. 

All the power spectral density plots (Figures 5.13 - 
5.18) exhibited a peak in power between 1x10° - 5x10-° Hz, 
corresponding to 1-2 days. This illustrates that most of 
the power is at the synoptic period. The typical slope of 
the psd function at the higher frequencies (above 10°° Hz) 
corresponding to periods less than one day was -1.5 in log- 


log space. 
E. ACCURACY OF HEAT FLUX MEASUREMENTS 


The largest source of uncertainty for SHF and LHF 
measurements 1S sensor accuracy (Blanc 1987). Using the 
manufacturers published sensor errors along with field 
experience (Table 5.9) an error estimate (Tables 5.10 and 
>.11) caused by sensor inaccuracies, ship contamination and 
radiation contamination for our exercises was determined in 


a root mean square manner (Equations 5.1 and 5.2). 
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Table 5.9 Possible measurement errors of data used in heat 








flux calculations. 


SHF,,, = (Pp CoC, (AT,,,) (U))#+ (p CoC. (Ad) (U3) )* (5.4) 


LHF.,,=y (p L,C,(Aq,,,) (U))?+ (pL,C,(Aq) (U,,,))?% (5.2) 
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Table 5.10 Possible error in sensible heat flux (W/m?) 


measurement as a function of wind speed and ASTD. 






Table 5.11 Possible error in latent heat flux (W/m?) 


measurement as a function of wind speed and air-sea humidity 


diference. 


The random SHF and LHF errors for average values of 
wind speed (10 m/s), ASTD (-5 C) and air-sea humidity 
difference (-1.5 g/Kg) from the NPS data set would be 21 
W/m? and 18 W/m’ respectively. Therefore the random THF 


DZ 








error would be about 27 W/m’, when added in a root mean 
Square fashion. Random errors can however be reduced by 
increasing the number of averaging operations done. This 
can be accomplished without changing method by taking more 
measurements. Systematic (non-random) errors are also 
present and cannot be reduced by averaging. The worst case 
non-random error for the NPS exercises would be 39 W/m. 
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Figure 5.1 Time series of sensible, latent and total 
turbulent heat flux (W/m) for SIZEX in the Greenland Sea, 
January 1992. 
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Figure 5.2 Histogram of total turbulent heat flux (W/m) for 
SIZEX in the Greenland Sea, January 1992. 
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Figure 5.3 Time series of sensible, latent and total 


turbulent heat flux (W/m*) for CEAREX in the Greenland Sea, 
March 1989. 
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Figure 5.4 Histogram total turbulent heat flux (W/m*) for 
CEAREX in the Greenland Sea, March 1989. 
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Figure 5.5 Time series of sensible, latent and total 
turbulent heat flux (W/m’) for NORCSEX88 in the Norwegian 


Sea, March 1988. 
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Figure 5.6 Histogram of total turbulent he 


NORCSEX in the Norwegian Sea, March 1988. 
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Figure 5.7 Time series of sensible, latent and total 
turbulent heat flux (W/m*) for NORCSEX in the Norwegian Sea, 


November 1991. 
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Figure 5.8 Histogram of total turbulent heat flux (W/m?) for 
NORCSEX in the Norwegian Sea, November 1991. 
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Figure 5.9 Time series of sensible, latent and total 


turbulent heat flux 
February-March 1989. 


(W/m) 


for CEAREX in the Barents Sea, 
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Figure 5.10 Histogram of total turbulent heat flux (W/m) 
for CEAREX in the Barents Sea, March 1989. 
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Figure 5.11 Time Series of sensible, latent ant total 


turbulent heat flux (W/m*) for SIZEX in the Barents Sea, 
March 1992. 
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Figure 5.12 Histogram of total turbulent heat flux (W/m*) 
for SIZEX in the Barents Sea, March 1992. 
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Figure 5.13 Power spectral density of latent (dashed LIne): 7 
sensible (dotted line) and total turbulent (solid line) heat 
Flux (W/m*) for SIZEX in the Greenland Sea, January 1992. 
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Figure 5.14 Power spectral density sensible (dotted line), 
latent (dashed line) and total turbulent (solid line) heat 
Flux (W/m*) for CEAREX in the Greenland Sea, March 1989. 
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Figure 5.15 Power spectral density of sensible (dotted 
line), latent (dashed line) and total turbulent (solid line) 
heat flux (W/m*) for NORCSEX88 in the Norwegian Sea, March 


Loses. 
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Figure 5.16 Power spectral density of sensible (dotted 
line), latent (dashed line) and total turbulent (solid line) 
heat flux (W/m’) for NORCSEX91 in the Norwegian Sea, 
November 1991. 
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Figure 5.17 Power spectral density of sensible (dotted 
line), latent (dashed line) and total turbulent (solid line) 
heat flux (W/m*) for CEAREX in the Barents Sea, February- 
March 1989. | 
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Figure 5.18 Power spectral density of sensible (dotted 
line), latent (dashed line) and total turbulent (solid line) 
heat flux (W/m*) for SIZEX in the Barents Sea, March 1992. 
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VI. CONCLUSIONS 


This thesis examined the exchange of heat between the 
ocean and atmosphere in the Greenland, Norwegian and Barents 
seas. The bulk method was used to calculate the turbulent 
heat fluxes in the winter based on bulk data gathered aboard 
ships. The results showed the largest positive fluxes 
occurring in the Greenland Sea, followed by the Norwegian 
sea and the smallest fluxes in the Barents Sea. The 
Greenland Sea is also a location of strong convective 
overturning and deep water formation in the Arctic. 

The values of the turbulent heat flux were generally 
found to be consistent with previous climatological results. 
However, HC flux estimates in proximity to the ice edge were 
much larger than those observed during the NPS cruises. 
Specifically, they estimated mean turbulent heat fluxes well 
in excess of 400 W/m’? and even as high as 600 to 700 W/m? in 
portions of the Greenland and Barents Seas. Based on the 
NPS data and in comparison to the climatological works in 
the regions, the estimates of HC are extremely high. In 
fact their mean values for some regions of the Greenland and 
Barents Seas are nearly equivalent to maximum fluxes found 
in this study. Differences in their heat flux values and 
those in this thesis may be due in part to larger winds in 
the HC data. Also, HC used moisture and heat exchange 
coefficients about 11% and 333 larger than our measured 
coefficients, explaining some of the differences. 

A large uncertainty in the flux measurements of Bunker 
and HC resulted from the use of questionable heat and 
moisture exchange coefficients. Bunker used a stability and 
wind dependent exchange coefficients that appeared to give 
values about 40% larger than for the method of Smith (1988). 


vee 





HC used constant exchange coefficients that would be 11-33% 
larger than those of Smith (1988) for the NPS data sets used 
in this thesis. 

Finally, the averaged values of the bulk meteorological 
parameters can be used in the bulk formula to obtain ~ 
averaged sensible and latent heat fluxes within 5% of the 
actual values. These results validate the classical monthly 
averaged method (SAM) employed by VT, Gorshkov and others in 
the Arctic Seas. However, other errors in their method may 
cause inaccurate results. A principal advantage of the 
classical method is speed of processing. After the raw data 
has been quality controlled, an average heat flux can be 
computed in seconds via the classical method compared with 


one hour or more with the sampling method. 
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APPENDIX. [MATLAB PROGRAM FOR COMPUTING HEAT FLUX] 


matlab 4.0 program run.m 
updated Feb 1995 
based on smith (jgr,v 93,n0 cl2, 15467-15472, 1988) 
Calculate the following fluxes 


shf sensible heat flux 
lht latent heat flux 
hf (total) heat flux 
based on | 
rho atm denstiy based on sfc temp 


ustar see below 
tstar see below 


N oo d dP oP oP dO oP OP GO dO oO oO oO ClO 


dy specific heat set constant at 1004 j kg-1 k-1 
Cc latent heat of vaporiz const at 2.50e+6 j kg-1 
CSSSSSSSSSESESESSESESSSESESSSSESYELLY LAL YYYIL AI 9 999999999009 
USl5us 'S 
ZS Su, 
load norx88; 
myfile=norx8s; 
[n m]=size(myfile) ; 
for |=iens 
Minty) = el) * 20s 
date(j) amyfile(j,1} ; 
time (j)=myfile(j,2); 
utrue(j)=myfile(j,7); 
tair(j)=myfile(j, a. 
tsfic(j)=myfile(j,4); 
hig) ae, 5)4 
p(j)=myfile(j,6); 
wdt (j) =myfile(j, 8) ; 
lat(j)=myfile(j,9); 
lon (3) =myfile (3.10) ; 
jd(j) = 67 + (1200 / 1440) + ( min(j) / 1440); 
wind(j)=utrue(j); 
astd(j)=tair(3)- =tste(7)}: 
press (j)=p(3); 
relhum(j)=rh(j); 
Sst) )=tsters). 
temp (j)=tair(3j); 
utrue=utrue(j); 
talr= tair(9); 
CSta= tsfic(j Ls 
rh=rh(4).; 
SoG) 
S calculate bulk ustar,tstar,qstar 
S$ anput: utrue wind speed (m/s) at zu 
% tair alr temperature (centigrade) at zt 


rh relative humidity (%) at zt 


ips 


HW AP AP AP DO AO AP oP GO GO dO DO GO oO CO 


oo 


HW dO oP oP oO oO 


re )\9) 





ESEC surface temperature (centigrade) 


ZU measurement level of utrue (centigrade) 
Ze measurement level of tair and rh (m) 
p atmospheric pressure at zt (mb) 


OUEDUC: 

ustar friction velocity (m/s) 

tstar scaling temperature (k or c) 

qstar scaling specific humidity (g/kg) not (g/g)! 

il monin-obukov length scale (m) 
calculates tstar based on a constant chn,cen 
zo base on charnocks relation 
this follows smith(1988) as close as psby except different 
virtual temperature calculation is used (smith was wrong) 
the surface 1s assumed to be saturated. 
set constants 


Z1LO0=10; ' %! reference height (not measuremement height) 
next two are based on value at 210 
ectn=1..0e=3- 6! heat and buoyancy flux xfer param at 210 
cen=1.2e-3; $! humidity flux transfer parameter 
K=.4> $! von karmen’s constant 
g=9.8; $! gravity 
a=0.011; 6! charnocks’s constant from smith(1980) 
b=0.1 1 6! smooth flow constant from businger (1973) 
nu=1.4e-5; 6! dynamic viscosity of air 
gamma=0.0096; %! adiabatic lapse rate 
SESCCESESESESESESEEEESESESSSSSEESESEESESEESSSSSSESESSSESSSYLLLLELES 
calculate parameters not requiring iteration 
( 


calculate saturation mixing ratio, qsat 
based on temperature (c) and press (mb) from stull (1988) 
it 1s only valid near the surface.use loew’s polynomial for 
more accurate or upper level estimates of qsat. 
psfc=p+0.116 * zt; % sfc press (mb) based on standard atms 
es=6.1078*exp (17.2694 *tsfc/ (tsfic+237.3)):% vap p(mb) tsfc 


qsat=622.0 * es/ (psfc-es) ; $ (g/kg) 
esa=6.1078*exp (17.2694 *tair/ (tair+237.3));%vap p(mb) tsfc 
qair=622.0 * esa/ (p-esa); 6 (g/kg) 


q=qair * rh/100; % mixing ratio at zt (g/kg) 

qsfc=qsat; %! mixing ratio at surface assumed staturated 
theta=tair+273.16+gamma * zt; %! potent temp (k) at zt 

the following is based on stull (1988) ;smith(1988) is wrong 
thetav=theta* (1.0+0.61e-3*q); %! virt potent temp (k) at zt 


tvsfc=(tsfic+273.16)*(1.0+0.61le-3*qsfc) ;% v 
SCSCCCCECECESESSESESEESESESESESESSESESELELEL SSS 
initialization (these values will change) 
count=0% 
fe). 0eEl0: 
ZOo=1 .0e-4; 


ustar=utrue * .036; 
find ustar, tstar loop; exit after convergence or 40 cycles 
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6 count is var for the iterations 
1f utrue > 0.2; 
calculate an init ustar 
ustarl=ustar; 
Zuko fF 7s 
ZU = zt. f As 
ratio = zu: / -zo-; 
calculate the integral diabatic heating term 
for momentum. we have used values from dyer (1974) 
unstable 
1f zu/l < 0.; 
Kelso 16.. *(ZU7 1) \ 225% 
psim=2.*logm((1.+x)/2.)+logm( (1.+x*x) /2)-2*atan(x)+3.14159/2: 
6 stable 


oo 


oP oO oO 


else 
psim = -5 * (zu/ 1); 
end 
* for temperature. based on large and pond (jpo v. 12,1982) 
6 unstable 


LE 2u/) < O44 
Kelle 16." (27 1) ) 052 
psit=2.*logm((1.+x*x) /2); 

stable 

else 
psit' = -5 * (zgu/ 1); 

end 

if ratio < 5.0 ; % prevent underflow 
count=40; 

end 

zOot=zl0*exp(-k*k / (ctn * logm (ZU7Z0}.))'s & 

Cstarv= (thetav-tvsfc) * k / (logm(zt / zot)-psit); 

ustar=utrue* k / (logm (zu / zo)-psim) ; 


oe 


zc=a * ustar* ustar/g; % charnocks formula 
ZS =b*nu / ustar; * smooth formula 
ZO=ZC+ZS; 
if abs(tstarv) < 1.0e-20; 

1l=1.0e10; 
else 

l=ustar* ustar* tvsfc/ (g*k * tstarv) ; 
end 


do iterations to find ustar convergence 
if abs((ustar-ustarl)/ustar) > .0005; 
for i=1:40; 
ustarl=ustar; 
20k. = 20..f Jos 
Zed = Ze. As 
ratio = Zu. / -zo% 
this function calculates the integral diabatic heating term 
for momentum. we have used values from dyer (1974) 


oo 


oe oO 


cae 


% 





unstable 
2 uy). =< .0e 
Sa Ll 6 (es) 


psim=2.*logm((1.+x)/2.)+logm( (1.4+x*x) /2)-2*atan(x)+3.14159/2; 


2 
xe) 


oP oO 


fo) ed 


stable 
else 
DSsim-= <=5.* (zu/ 2) 
end 
for temperature. based on large and pond (jpo v. 12,1982) 
unstable © 
1f 2u/l = 0. 
x] Cl HO. (2071) 225 
psit=2 <* logm((1.4x*x)/2)% 
stable 
else 
psit = =5 * (2a 2 
end 
if ratio < 5.0 ; % prevent underflow 
count=40; 
end 
zot=z10*exp(-k*k / (ctnm * Logm (zu/zo))); & 
tstarv= (thetav-tvsfc) * k / (logm(zt / zgot)- 
ustar=utrue* k / (logm (zu / zo)-psim) ; 
zc=a * ustar* ustar/g; % charnocks formula 
ZS =b*nu / ustar; * smooth formula 
ZO=ZC+ZS; 
if abs(tstarv) < 1.0e-20; 
1l=1.0e10; 
else | 
l=ustar* ustar* tvsfc/ (g*k * tstarv); 
end 
end %$ is for iat 40 
end %@ for abs ( ) > 
end % when u<.2 
post loop calculations 
if count < 40 & abs((ustar-ustarl)/ustar) < .0005; 
tstar=(tair-tsfc) * k / (logm(zt/ zot)-psit); 
qstar=(q-qsfc) *k / (logm(zt/ zot)-psit); 
else @ too stable loop does not converge 
uUStar=U..2 
tstar=0.; 
Gstare0..; 
t=O) .% 
end 
rho=psfc/(2.87*tvsfc); % calc var density with ideal gas law 
shf(j)=-rho*1004 * ustar* tstar; 
lnf (3) =-rho*2.5e6 * ustar* qstar/1000;1000 
HE (5) =sht(3) +1 ht 


pSsit); 
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cd(j)=(uustar(j)/ utrue) * 2; 

cdn(j)= (1 / sqrt(cd(j))+psim/.4) * (-2); 

fprintf (‘working on observation number %g\n’,j) 
end % stop counting j here 
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